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ABSTRACT 

Using a near-IR stellar library of 706 stars with a wide coverage of atmospheric 
parameters, we study the behaviour of the Can triplet strength in terms of effective 
temperature, surface gravity and metallicity. Empirical fitting functions for recently 
defined line-strength indices, namely CaT*, CaT and PaT, are provided. These func¬ 
tions can be easily implemented into stellar populations models to provide accurate 
predictions for integrated Can strengths. We also present a thorough study of the 
various error sources and their relation to the residuals of the derived fitting func¬ 
tions. Finally, the derived functional forms and the behaviour of the predicted Can 
are compared with those of previous works in the field. 

Key words: stars: abundances - stars: fundamental parameters - globular clusters: 
general - galaxies: stellar content. 


1 INTRODUCTION 

This is the third paper in a series dedicated to the under¬ 
standing of stellar populations of early-type galaxies and 
other stellar systems by using their near-IR spectra and, 
particularly, the strength of the integrated Call triplet. In 
the previous papers, we have presented the basic ingredi¬ 
ents of this project. First, we have observed a new stellar 
library of 706 stars at 1.5 A (FWHM) spectral resolution in 
the range AA 8348-9020 A (Cenarro et al. 2001a, Paper I). 
That paper also includes the definition of 3 new, improved 
line-strength indices in the Ca ll triplet region (namely CaT, 
PaT, CaT*), which are especially suited to be measured in 
the integrated spectra of stellar populations. Paper II of the 
series (Cenarro et al. 2001b) presents an updated set of at¬ 
mospheric parameters (Teff, logg and [Fe/H]) for the stars 
of the library. The objective of this third paper is to pro¬ 
vide empirical fitting functions describing the behaviour of 
the above indices in terms of the atmospheric parameters. 
Together with the spectra of the stellar library, the fitting 
functions will be implemented into an evolutionary stellar 
populations synthesis code to predict both the spectral en¬ 
ergy distribution and the strengths of the indices for stellar 
populations of several ages and metallicities (Vazdekis et al. 
2001, Paper IV). 

* E-mail: cen@astrax.fis.ucm.es 


A popular method to investigate the star formation his¬ 
tory and the stellar content of galaxies is to compare ob¬ 
served line strength indices with stellar population models. 
This requires an accurate, prior knowledge of the behaviour 
of the spectral features of interest for a wide range of stellar 
spectral types, luminosity classes and metallicities, which 
can be accomplished with the aid of either empirical stellar 
libraries or theoretical model atmospheres. There are two 
common alternatives for the stellar populations synthesis 
approach. The first method mixes the spectra of the differ¬ 
ent stars with their relative ratios given by evolutionary syn¬ 
thesis models to predict spectral energy distributions (Fioc 
& Rocca-Volmerange 1997; Leitherer et al. 1999; Vazdekis 
1999; Schiavon, Barbuy & Bruzual 2000; Bruzual & Chariot 
2001). The major advantage of this spectral synthesis is that 
it provides full information of all the spectral features within 
the spectral range covered by the stellar library. However, 
its usefulness relies on the availability of a complete stellar 
library at the appropriate spectral resolution. The second 
procedure, and probably the most widely employed in the 
past, is the use of empirical fitting functions describing the 
strength of previously defined spectral features in terms of 
the main atmospheric parameters. These calibrations are di¬ 
rectly implemented into the stellar populations models to 
derive the index values for populations of different ages and 
metallicities. It must be noted that one of the main advan¬ 
tages of using fitting functions to reproduce the behaviour 
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of spectral indices is that it allows stellar populations mod¬ 
els to include the contribution of all the required stars by 
means of smooth interpolations between well-populated re¬ 
gions in the parameter space. One should not forget that, 
in any of the above two methods, a common limitation of 
the empirical procedures arises from the fact that they im¬ 
plicitly include the chemical enrichment history of the solar 
neighborhood. This caveat must be kept in mind when us¬ 
ing model predictions to interpret the stellar populations of 
external galaxies, whose star formation histories might be 
completely different to that of the Galaxy. 

The usefulness of the fitting functions approach has 
been clearly demonstrated by the fact that the most im¬ 
portant evolutionary synthesis models (e.g., Worthey 1994, 
Vazdekis et al. 1996, Tantalo et al. 1996 or Bruzual & Char¬ 
iot 2001) have implemented the available fitting functions to 
fit observed line strengths in the literature. Unfortunately, 
at present fitting functions are only available in the blue and 
visible part of the spectrum (e.g. Gorgas et al. 1993, Worthey 
et al. 1994 and Worthey & Ottaviani 1997 for the Lick/IDS 
indices; Poggianti & Barbaro 1997 and Gorgas et al. 1999 
for the A 4000A break). Since other spectral regions also 
contain very useful, complementary, absorption lines, it is 
necessary to extend this kind of calibrations to line-strength 
indices in other spectral regimes, such as the ultraviolet and 
the near infrared (e.g. the CO index at 2.2 pm analyzed by 
Doyon, Joseph & Wright 1994). This is indeed an important 
motivation for this paper. So, before we give a comprehen¬ 
sive analysis of the predictions of both the spectral synthesis 
and the fitting functions for stellar populations (Paper IV), 
we will concentrate in this paper on understanding the be¬ 
haviour of the Call triplet as a function of the atmospheric 
parameters, providing the corresponding fitting functions. 

The Call triplet is one of the most prominent features 
in the near-IR region of the spectrum of cool stars. Even 
though there are many previous works dealing with the near- 
IR Ca II triplet (we refer the reader to the short review pre¬ 
sented in Section 2 of Paper I), the papers by Dfaz, Ter- 
levich & Terlevich (1989), Zhou (1991), Jprgensen, Carlsson 
& Johnson (1992) and Idiart, Thevenin & de Freitas-Pacheco 
(1997) (hereafter DTT, ZHO, JCJ and ITD, respectively) 
deserve a special mention, since they provide empirical fit¬ 
ting functions for the Call triplet. They have been very 
useful to obtain a first understanding of the behaviour of 
the Call triplet as a function of the stellar parameters but, 
as we will see in this paper, suffer from several limitations 
that make it impossible for stellar populations models to 
predict reliable calcium strengths, especially for old stel¬ 
lar populations. Previous papers which have made use of 
the above mentioned fitting functions to predict the Call 
triplet in the integrated spectra of galaxies include Vazdekis 
et al. (1996), ITD, Mayya (1997), Garcfa-Vargas, Molla & 
Bressan (1998), Leitherer et al. (1999) and Molla & Garcfa- 
Vargas (2000). 

Section 2 of this paper describes the qualitative be¬ 
haviour of the Call triplet as a function of the atmospheric 
parameters, as derived from the new stellar library. We also 
include a comparison with the results of the previous work. 
We devote Section 3 to the mathematical fitting procedure, 
providing the significant terms, coefficients and statistics of 
the derived fitting functions. Afterwards, a thorough anal¬ 
ysis of residuals and possible error sources is presented, in- 


Table 1. Bandpass limits for the generic indices CaT and PaT. 


CaT central 
bandpasses (A) 

PaT central 
bandpasses (A) 

Continuum 
bandpasses (A) 

Cal 8484.0-8513.0 
Ca2 8522.0-8562.0 
Ca3 8642.0-8682.0 

Pal 8461.0-8474.0 

Pa2 8577.0-8619.0 
Pa3 8730.0-8772.0 

8474.0-8484.0 

8563.0-8577.0 

8619.0-8642.0 

8700.0-8725.0 

8776.0-8792.0 


eluding the sensitivity of the fitting functions to differences 
in Ca/Fe ratios. In Section 4, we compare the new fitting 
functions with those presented in previous papers. Finally, 
Sections 5 is reserved to discuss some important issues and 
to summarize the contents of this paper. 


2 BEHAVIOUR OF THE CA II TRIPLET AS A 

FUNCTION OF ATMOSPHERIC 

PARAMETERS 

2.1 Qualitative behaviour 

As a first step to understand the behaviour of the Call 
triplet as a function of the atmospheric parameters, this 
section describes, from a qualitative point of view, the ef¬ 
fect of temperature, surface gravity and metallicity on the 
strength of the Call lines. Figure ^ shows a comparative 
sequence in spectral types for a sample of dwarfs (a) and 
giants (b) taken from the new near-IR stellar library (see 
also Figure 1 in Paper I for a detailed description of the 
strongest spectral features). The effects of metallicity and 
gravity on the Call lines are clearly shown in Figure ^ In 
Figure ^ we plot the measured CaT, PaT and CaT* indices 
versus 9 (= SOdO/Teii) for the whole sample of stars. 

It is worth reminding that the indices CaT and PaT 
measure, respectively, the strengths of the raw calcium 
triplet and of three pure H Paschen lines. They consist of 
five continuum bandpasses (which are the same for both in¬ 
dices) and three central bandpasses covering, in each case, 
the three Call lines (CaT) or the P17, P14 and P12 lines 
of the Paschen series (PaT). Table ^ lists the bandpass 
limits of these indices. The PaT index was specially de¬ 
signed to quantify the CaT contamination with Paschen 
lines, and it makes it possible to define the corrected index 
CaT*(= CaT — 0.93PaT) as a reliable indicator of the pure 
Call strength. We refer the reader to Paper I for further 
details about the definition and measurement of the above 
indices. 

The spectra of hot stars (Fig. ^ top) are dominated by 
Paschen lines, whose relative depths exhibit a smooth se¬ 
quence with wavelength, with the Call strengths, both in 
dwarfs and giants, being clearly negligible (see also Fig. |c). 
As temperature decreases, the Can lines become prominent, 
standing out over the Paschen lines sequence (see the A and 
F types in Fig. The increasing strength of the Call index 
with the decreasing temperature peaks at late K and early 
M types. For a wide range of spectral types (0.7 < 9 < 1.2), 
the Ca II lines are heavily affected by metallicity and gravity 
effects (leading to the large spread of CaT and CaT* values 
in Figs. ^ and ic) , in the sense that their strengths increase 
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wavelength (A) wavelength (A) 

Figure 1. Sequences in spectral types for (a) dwarf and (b) giant stars from the near-IR stellar library. Effective temperatures, names, 
spectral types, luminosity classes and surface gravities (log g in dex) are given in the labels. All the spectra have been normalized and 
reproduced using the same scales so relative differences among the spectra are kept. 


as metallicity increases and gravity decreases (Fig. |^. For 
lower temperatures, Ca goes gradually into a neutral state 
and the Call strength decreases. The turning point at which 
the trend with temperature changes depends on the lumi¬ 
nosity class, being cooler for giants (9 ~ 1.3) than for dwarfs 
{9 ~ 1.1). Such an effect is readily seen in Figs. ^ and 
and also by comparing the depths of the Call lines for the 
M4 V and M4 III types in Fig. |^. Finally, strong molecular 
bands of TiO and VO dominate the latest spectral types, 
in which the Call strength is almost negligible (Fig. bot¬ 
tom). For these stars, as well as for the very hot ones, metal¬ 
licity seems not to affect the calcium strength. 

The behaviour of the PaT index is illustrated in Fig ^b. 
As expected for a hydrogen index, it attains high values for 
hot and warm stars {9 ^ 0.8), showing a maximum at an 
effective temperature which depends on the luminosity class 
(already reported by Andrillat, Jaschek & Jaschek 1995). 
Even though, in the low temperature regime, the stellar 
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spectra do not present Paschen lines, the index takes values 
slightly larger than zero. This is mainly due to the contami¬ 
nation of weak metallic lines falling within the central bands. 
Also, the presence of strong molecular bands in the spectra 
of cold giants {9 ^ 1.4) causes a fictitious enlargement of 
the PaT values. 


2.2 Comparison with previous work 

In Paper I we have already presented a brief review of the 
previous work studying the Ca ll triplet as a function of at¬ 
mospheric parameters. To facilitate the following discussion, 
in Table ^ we give a summary of the conclusions of these 
previous studies. We include the main empirical papers on 
the subject, together with the theoretical calibration of JCJ. 
Columns 3 to 5 give the ranges of the atmospheric parame¬ 
ters spanned by the calibrating stars (the numbers of stars 
used are listed in column 2) or atmosphere models. The at- 
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wavelength (A) wavelength (A) 

Figure 2. Metallicity and gravity effects on the strength of the Can triplet for a subsample of stars from the near-IR stellar library. 
Panel (a) shows stars with similar temperature and gravity but spanning a wide range in metallicity. Panel (b) displays a sequence in 
gravity for stars with similar temperature, and metallicity around solar. Temperatures in K are given in brackets. All the spectra have 
been normalized and reproduced using the same scale. 


mospheric parameters on which the Call triplet was found 
to be dependent in these previous works are indicated in 
columns 6 to 8, where the keywords strong and weak are 
used to remark that the Call strength was reported to be 
more sensitive to one parameter than to another. 

A first glance at Table ^ reveals that there are several 
apparent inconsistencies among the previous papers. These 
discordant results are mainly due to differences in the cov¬ 
erage of the parameter space and the index definitions. One 
of the points of discrepancy is the sensitivity to effective 
temperature. For the low temperatures, as it was first re¬ 
ported by Cohen (1978), Figs. and ^ clearly show that 
the Call strength decreases for stars colder than 9 ~ 1.3. 
This behaviour has been thus acknowledged by the papers 
which include a good number of M stars, like Carter, Vis- 
vanathan & Pickles (1986) and ZHO. However, other au¬ 
thors, in particular Jones, Alloin & Jones (1984; hereafter 


JAJ) and DTT, were not able to detect such an effect. The 
problem here is that the index definitions that they used are 
not well suited for stars with strong TiO absorption bands 
(see both Section 4.2 and Figure 2 of Paper I). In the first 
case (JAJ), the index is uncertain for cold stars since the red 
sideband falls in a strong TiO absorption. In the second one 
(DTT), the position of the local continuum causes the index 
to overestimate the Ca ll absorption. Concerning the hotter 
stars. Fig. ^ shows that the Call strengths decrease when 
Teff increases. However, and with the exception of Alloin & 
Bica (1989; hereafter A&B), this behaviour has not been 
noted in the previous empirical papers since the Call lines 
are heavily contaminated by the Paschen series for 6 < 0.8 
(compare Figs. and |c). Note that JCJ do reproduce this 
dependence since they compute synthetic calcium equivalent 
widths without the inclusion of Paschen lines. 

Concerning the gravity dependence. Table ^ shows that 
© 2000 RAS, MNRAS 000, 0-0 
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Figure 3. CaT, PaT and CaT*(= CaT — 0.93PaT) indices versus 6 {= 5040/Teff) for the whole stellar library. Different symbols are 
used to indicate different metallicities (as in the key), while sizes are related with surface gravity, in the sense that the larger the symbol 
(supergiants), the lower the gravity. On the top, the effective temperature scale is given. 
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Table 2. Previous calibrations of the Call triplet. Sources are the following: JAJ (Jones, Alloin & Jones 1984), CVP (Carter, Visvanathan 
&; Pickles 1986), A&B (Alloin & Bica 1989), DTT (Diaz, Terlevich &; Terlevich 1989), ZHO (Zhou 1991), MAL (Mallik 1997), ITD (Idiart, 
Thevenin &; de Freitas Pacheco 1997), and JCJ (Jdrgensen, Carlsson & Johnson 1992). 


Source 

No. 

stars 

e 

range 

logff 

range 

[Fe/H] 

range 

Teff y [Fe/H] 

dependences 

Fitting functions 
parameters 

JAJ 

62 

0.68 - 

1.47 

0.7 - 4.8 

-0.60 - 

-1-0.43 

no 

strong 

weak 


CVP 

51 

0.93 - 

1.68 

0.0 - 4.6 

0 - 

-1-0.50 

yes 

yes 

no 


A&B 

63 

0.68 - 

1.46 

0.7 - 4.8 

-0.60 - 

-1-0.43 

? 

yes 

yes 


DTT 

106 

0.74- 

1.43 

0.2 - 4.6 

-2.70 - 

-1-0.55 

no 

yes 

yes 

logy, [Fe/H] 

ZHO 

144 

0.76 - 

1.72 

0.7 - 4.8 

-2.28 - 

-1-0.60 

yes 

yes 

yes 

(R-/), {R-If, logy, log^y, lO^'^/Hl 

MAL 

146 

0.81 - 

1.58 

-0.6 - 4.8 

-3.0 - 

-1-1.01 

weak 

yes 

yes 


ITD 

67 

0.74- 

1.20 

0.8 - 4.5 

-3.15 - 

-1-0.35 

yes 

weak 

strong 

0, logy, [Fe/H], 0[Fe/H], logy[Fe/H], [Ca/Fe] 

JCJ 

- 

0.76 - 

1.26 

0.0 - 4.0 

-1.0 - 

-1-0.20 

yes 

yes 

yes 

Feff, Tff, logy, log^y, T^b logy 


most previous papers have remarked the strong luminosity 
sensitivity of the Call triplet. One exception is the work by 
ITD, in which they found a weaker dependence on gravity. 
However, it should be noticed that their sample comprises 
only a small number (8) of supergiants, and that most of 
them (7) exhibit very low metallicities ( — 2.6 < [Fe/H] < 
— 1.4). In other words, gravity and metallicity effects may 
not be well separated. Note also that, as it is shown be¬ 
low, gravity effects at the low metallicity range are milder. 
The discrepancies among different works, as far as the sen¬ 
sitivity to metallicity is concerned, are clearly due to the 
[Fe/H] range spanned by the calibrating stars. While the 
early works by JAJ and Carter et al. (1986) could not de¬ 
tect a significant metallicity dependence, subsequent papers 
which explored a broader [Fe/H] range (like DTT and ITD) 
clearly showed that metal abundance is also a key parame¬ 
ter. 

2.3 The gravity and metallicity dependences 

In order to explore in more detail the interrelationship be¬ 
tween gravity and metallicity effects on the Call triplet of 
the cool stars, we plot in Figure ^ the gravity dependence 
for different metallicity ranges, for stars with effective tem¬ 
peratures 0.8 < 9 < 1.3. To guide the eye, we have fitted 
parabola’s for each metallicity bin (note that these lines are 
not the definitive fitting functions), checking that the sec¬ 
ond order term is statistically signihcant in all cases. From 
this figure, it is clear than the dependence on gravity is not 
linear. This shows that the Call triplet is much more sen¬ 
sitive to gravity for giants and supergiants than for dwarfs 
(for logy ^ 3.5 the relation is almost flat). This behaviour 
contrasts with the results of JAJ, A&B and DTT, who find a 
linear relation with logy, but it is in qualitative agreement 
with ZHO and, especially, with the conclusions of Mallik 
(1997; hereafter MAL). Another result of Fig. ^ is that the 
gravity dependence tends to flatten off for low metallicities. 
This effect was already noted by DTT when they concluded 
that, at low Z, metallicity was the main parameter and it 
is, again, in perfect agreement with the results of MAL (see 
their Figure 6 and, also, Mallik 1994). It must also be noted 
that our conclusions about the gravity dependence also agree 
extremely well with the qualitative predictions of the the¬ 
oretical study by JCJ. We even notice their result about 
a Call increase at high gravities for metallicities around 
[Fe/H] = -1. 



Figure 4. Gravity dependence of the CaT* index for stars in the 
temperature range 0.8 < 0 < 1.3. Symbols types indicate different 
metallicities following the code of Fig. |^. The lines represent 
quadratic fits in logy for stars in the metallicity ranges: [Fe/H] > 
0.25, -0.25 < [Fe/H] < 0.25, -0.75 < [Fe/H] < -0.25, -1.25 < 
[Fe/H] < -0.75, -1.75 < [Fe/H] < -1.25, and [Fe/H] < -1.75 
(from top to bottom). 


The metallicity dependence is analysed in more detail 
in Figure ^ where we have grouped our cool stars into three 
luminosity classes. We also show linear fits to the [Fe/H] 
dependence of each gravity subsample (in this case, the sec¬ 
ond order terms are not statistically significant). The CaT*- 
metallicity relation is thus roughly linear for the cool stars, 
in agreement with A&B. Furthermore, we must note that we 
do not find any flattening of the metallicity relation at high 
[Fe/H], as was reported by DTT (although they derived a 
linear relation with [Fe/H], they concluded that, in the high 
metallicity range, the Call strength is a function of grav¬ 
ity only). More interestingly, we Hnd that the slope of the 
metallicity relation is statistically the same for giants and 
dwarfs (already noticed by A&B), but much steeper for the 
supergiant subsample. Again, this is in agreement with the 
findings of MAL. In particular, at low metallicities, giants 
and supergiants do not differ in their Call strengths (as it 
was already apparent in Fig 

To summarize, we Hnd that the Ca ll triplet in cool stars 
(~ F7 - MO spectral types) follows a complex dependence 
on temperature, metallicity and gravity. Our conclusions are 
in very good agreement with the empirical results of MAL 
and the theoretical predictions of JCJ. 
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[Fe/H] 

Figure 5. Metallicity dependence of the CaT* index for stars 
in the temperature range 0.8 < 9 < 1.3. Different symbols are 
used to indicate different luminosity classes, as shown in the key. 
The straight lines are linear fits to each set of stars, with slopes 
of 3.6 ± 0.3, 2.2 it 0.1 and 2.0 it 0.1, for supergiants, giants and 
dwarfs, respectively 


3 THE FITTING FUNCTIONS 

In this section we present the empirical fitting functions for 
the new indices, as well as the general procedure followed 
to compute them. These fitting functions have been calcu¬ 
lated using the index measurements given in Paper I and the 
atmospheric parameters derived in Paper II. The complete 
tables of indices and parameters are available at the URL 


We chose the one that minimizes the residuals of the fit. 
Xa refers to any of the above indices and p is a polynomial 
with terms of up to the third order, including all possible 
cross-terms among the parameters, i.e., 


p(&dozg, [Fe/H]) = 


E 

0<i+j+k<3 


Ci.j.fcS* (log fif)-’[Fe/H] 


(3) 


with 0<i-|-i-l-fc<3 and 0 < i, j, k. 

Note however that, as a result of the wide parameter 
space covered by the stellar library, there exists no single 
function of the forms given by Eqs. (|^) and (^) which can 
accurately reproduce the complex behaviour of the indices 
CaT* and PaT. We have therefore divided the whole param¬ 
eter space into several ranges of parameters (boxes) in which 
local fitting functions can be properly computed. Finally, a 
general fitting function for the whole parameter space has 
been constructed by interpolating the derived local func¬ 
tions. To do so, we have defined the boundaries of the boxes 
in such a way that they overlapped, including several stars in 
common. In the overlapping zones cosine-weighted means of 
the functions corresponding to both boxes were performed. 
This would mean that, if Xai(x, y, z) and Xa, 2 {x, y, z) are two 
local fitting functions overlapping in the generic parameter 
z, and defined respectively in the intervals ( 21 , 1 , 21 . 2 ) and 
( 22 , 1 , 22 , 2 ) with 22,1 < 21 , 2 , the predicted index in the inter¬ 
mediate region will be given by 

l^{x,y,z) = wl^i{x,y,z) -I- (1 - w)I^ 2 {x,y, z), (4) 

where the weight w is modulated by the distance to the 


aaar esses: 

http://WWW.ucm.es/info/Astrof/ellipt/CATRIPLET.html 

overlapp 

ng limits as 

and 


TT / 2 - 22,1 \ 

http://www.nottingham.ac.uk/~ppzrfp/CATRIPLET.html 


2 1 ^ 21,2 - 22,1 J 


with 22,1 < 2 < 21 , 2 . 


(5) 


3.1 The fitting procedure 

The main objective of this paper is to derive empirical fitting 
functions for the indices CaT*, CaT and PaT in terms of the 
stellar atmospheric parameters. Keeping in mind that any of 
the above indices can be expressed as a linear combination 
of the other two ones, we have only computed fitting func¬ 
tions for CaT* and PaT. We chose them since these indices 
are better measurements of pure Call and H line-strengths 
respectively than CaT and, therefore, their behaviours with 
atmospheric parameters are expected to be simpler than in 
the case of a mixed line-strength like CaT. In any case, it 
is important to clarify that the empirical calibrations are 
only mathematical representations of the behaviour of the 
indices as a function of the atmospheric parameters and, 
thus, a physical justification of the derived coefficients is 
beyond the scope of this paper. 

Given the wide range of spectral types in the stellar 
library we prefer to use 6 as the effective temperature indi¬ 
cator, together with logg and [Fe/H], the classical parame¬ 
ters for surface gravity and metallicity. Following Gorgas et 
al. (1999) (see also Gorgas et al. 1993, and Worthey et al. 
1994), the fitting functions have been computed as polyno¬ 
mials of the atmospheric parameters in two possible func¬ 
tional forms, 

Xa(6l, logff, ]Fe/H]) =p(6', logg, [Fe/H]), or (1) 

T^iS,logg, [Fe/H]) = const. -I-exp [p{e,logg, [Fe/H])], (2) 


After trying several functional forms for the weights, we con¬ 
firmed that the analytical expression given in Eq. (^) guar¬ 
antees, in most cases, a smooth interpolation between local 
functions and preserves the quality of the final fit. 

The local fitting functions were derived through a 
weighted least squares fit to all the stars within each pa¬ 
rameter box, with weights according to the uncertainties of 
the indices for each individual star, as computed in Sec¬ 
tion 5 of Paper I. Since not all possible terms of Eq. (^ 
were necessary, we followed a systematic procedure to ob¬ 
tain the appropriate local fitting function in each case. To 
begin with, a general fit to all the 20 possible terms is com¬ 
puted together with the residual variance of the fit and the 
coefficients errors. The significance of each term is then cal¬ 
culated by means of a t-test (that is, using the error in that 
coefficient to check whether it is significantly different from 
zero) and the term with the highest significance level (q) is 
removed from the fit. The whole procedure in then iterated, 
removing in each turn the least significant coefficient, until 
all remaining coefficients are statistically significant. Typi¬ 
cally, we have used a threshold value of a = 0.10 to keep 
a significant term. It must be noted that the above proce¬ 
dure does not guarantee that we are obtaining the best solu¬ 
tion, so alternative approaches, like an inverse one (starting 
from a one parameter fit and introducing the most signifi¬ 
cant term in each iteration) were also followed to ensure that 
the final combination of terms was the one which provided 
the minimum unbiased residual variance. 

It should be noted that, after computing any local fit 
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Table 3. Stars which were not used for the fitting functions com¬ 
putation, coded as: C, carbon star; EmL, emission lines of the 
elements in brackets; P, pulsating star; SB, spectroscopic binary; 
Var, variable star; *: bad quality spectrum because of low signal- 
to-noise ratios, bad exposing conditions, unreliable spectral fea¬ 
tures, and others. 


Name 

Diagnostic 

Name 

Diagnostic 

HD 108 

EmL (Ca,H) 

HD 120933 

Var (CVn) 

HD 1326B 

Flare star 

HD 121447 

Var 

HD 17491 

p 

HD 138279 

Var (RR Lyr) 

HD 35601 

p 

HD 181615 

EmL (Ca) 

HD 42475 

p 

HD 217476 

Var 

HD 46687 

c 

HD 222107 

Var (RS CVn) 

HD 54300 

c 

BD+ 61 154 

EmL (Ca,H) 

HD 58972 

SB 

M5 IV-87 

* 

HD 60522 

Var 

M92 1-13 

* 

HD 74000 

* 

M92 H-23 

* 

HD 112014 

SB 

M92 VI-74 

* 

HD 115604 

Var 




in the above procedure, we made sure that the residuals did 
not present any systematic deviation, especially for stars of 
a given cluster or metallicity range. In some cases we found 
single stars with small index errors deviating excessively 
from the local fit, raising the residual variance of the final fit. 
When this occurred, the stars were analyzed in detail (look¬ 
ing for uncertainties in the atmospheric parameters, variabil¬ 
ity, chromospheric activity, spectroscopic binarity, etc) and, 
when necessary, thw were rejected from the fit. These stars 
are listed in Table M. Furthermore, all the previous steps in 
the fitting procedure (e.g. boxes design, overlapping zones 
definition, etc.) were also optimized to enlarge the quality 
and accuracy of the final fitting functions. 

3.2 Fitting functions for the indices CaT* and 
PaT 

The derived local fitting functions for the indices CaT* and 
PaT are presented, respectively, in Tables ^ and The Ta¬ 
bles are subdivided according to the atmospheric parame¬ 
ters ranges of each fitting box and include: the functional 
forms of the fits (polynomial or exponential), the signifi¬ 
cant coefficients and their corresponding formal errors, the 
typical index error for the N stars employed in each inter¬ 
val (utyp = N/ the unbiased residual variance 

of the fit ((Tgtd) 8.nd the determination coefficient (r^). Note 
that this coefficient provides the fraction of the index vari¬ 
ation in the sample which is explained by the derived fit¬ 
ting functions. Figure ^ displays the parameter regions (la¬ 
belled according to the atmospheric parameters regimes of 
Tables ^ and |) in which the final fitting function was calcu¬ 
lated by interpolating between the overlapping zones. Note 
that, in some cases, these boxes are somewhat narrower that 
the fitting regions listed in Tables and ^ This guarantees 
smoother transitions between different regimes. 

Users interested in implementing these fitting functions 
into their population synthesis codes can make use of the 
FORTRAN subroutine included in the URL addresses given 
above. This program performs the required interpolations 
to provide the CaT*, PaT and CaT (computed as CaT = 
CaT* -I- 0.93PaT) indices as a function of the three input 


30000 13000 7000 5000 4000 3350 3750 



Figure 6 . log 3 - 9 diagrams for the whole stellar library. The 
upper scale gives effective temperature in K. Different symbols 
are used to indicate stars of different metallicitie^ as coded in 
Fig. 1^. The boxes (labelled as in Tables ^ and display the 
regions of the corresponding local fitting functions for the CaT* 
(upper) and PaT (lower) indices. In the overlapping zones, the fi¬ 
nal fitting functions were derived by interpolating between boxes, 
as described in the text. Except for the box (c’) in the upper di¬ 
agram, specially designed to reproduce the proper behaviour of 
warm giants with low metallicities (dashed line), the rest of the 
regions are defined for the whole range of metallicities. 

atmospheric parameters. It also gives an estimation of the 
errors in the predicted indices (see below). 

3.2.1 CaT fitting functions 

As it is shown in Fig. the CaT* index in hot stars 
(0.13 0 0.69) shows a dichotomic trend for dwarfs on 

one side, and giants and supergiants on the other one. As a 
result of this, we designed boxes a and b in Fig. ^ (top) and 
derived different local fitting functions for the two groups 
of stars. Gravity effects are nearly negligible within each 
of the two luminosity bins. Also, due to the expected in¬ 
dependence on metallicity for high temperature stars, only 
terms in 9 were found to be statistically significant to re¬ 
produce the increasing CaT* trend with decreasing temper¬ 
ature (Tables ^ and^b). In the low temperature regime, 
the behaviour of cold dwarfs also differs from that of cold 
giants (Fig. ic). Here, we followed the same strategy and 
derived two different local fits in which, again, only terms 
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Figure 7. CaT* values and general fitting functions for different atmospheric parameter regimes. Panels a, b and c display, respectively, all 
the stars with gravities in the ranges 0.0 < logg < 2.0, 2.0 < logg < 3.2 and 3.2 < logg < 5.2, together with the derived fitting functions 
for the mean gravity in each range, i.e., logg = 1.0 (a), 2.6 (b) and 4.2 (c). In the mid-temperature range, the different lines represent, 
from top to bottom, the functions for metallicities [Fe/H] = -1-0.5, 0.0, -0.5, —1.0, -1.5 and -2.0, whilst, for high and low temperatures, 
the fitting functions do not depend on metallicity. Panel d shows all the stars around solar metallicity (—0.25 < [Fe/H] < -1-0.25) and 
the corresponding fitting functions computed for [Fe/H] = 0.0 and different values of gravity (logg = 0.0, 1.0, 2.0, 3.0, 4.0 and 5.0, from 
the thickest to the thinnest line). Codes and relative sizes of the star symbols (indicating, respectively, metallicity and gravity ranges) 
are explained in Fig. |^. Note that, while the lines displayed here correspond to fitting functions at particular values of logg and [Fe/H], 
the plotted stars span a range of atmospheric parameters around these central values. This is the reason why the lines do not exactly fit 
all the points in the plots. Also, note that these fitting functions have not been derived by using only these plotted stars, but the whole 
sample. 


in 6 were needed (Tables ^ and ^). As we mentioned in 
Section 2.1, the decreasing CaT* trend with increasing 9 


starts at lower temperatures for giants {9 ~ 1.30) than for 
dwarfs (9 ~ 1.07), being also steeper for the former group 
of stars. Finally, it is worth noting that, although a consid¬ 
erable number of hot and cold stars have unknown [Fe/H] 
determinations, the small scatter in the data shows that 
metallicity is indeed not governing the observed CaT* be¬ 
haviour at these temperature regimes. 

The complex behaviour of warm and cool stars was 
parametrized in terms of the three atmospheric parameters 
(see Tables ||;,|d and ^). Although an unique parameter 
box enclosing the rest of stars was enough to reproduce the 
CaT* behaviour at intermediate temperatures (box d in up¬ 
per panel of Fig. ^), we designed transition boxes (c and e) at 
both sides to improve the quality of the interpolations with 
the local functions derived for hot and cold giants stars. Note 
that cold dwarfs follow a trend smoother than cold giants. 
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and we did not need to include them in box d. Figures |^, 
^ and ^ illustrate the general fitting functions for three 
gravity bins (representing, roughly, supergiants, giants and 
dwarfs), showing an increasing dependence on metallicity 
as gravity decreases (therefore, supergiants span the widest 
range in CaT* values). Figure ^ displays the derived fitting 
functions for solar metallicity and different gravity values. It 
is clear from this figure that gravity effects increase with de¬ 
creasing gravity. Also, the higher the metallicity, the larger 
the gravity dependence. These effects were already apparent 
in Figs. ^ and |[ 

Due to the lack of warm metal-poor giants (there are 
only a few stars in the sample), the local fit in Table was 
specially designed to avoid unreliable interpolations between 
hot and warm giants in the low metallicity interval. More¬ 
over, in order to ensure a smooth interpolation, the limits 
of this box are slightly different depending on the range of 
metallicity (0.65 < 9 < 0.95 for [Fe/H] < —1.25, 0.65 < 9 < 
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Table 4. Coefficients and statistical data of the local fitting func¬ 
tions for the index CaT* in each range of atmospheric parameters. 
The term “giant” also includes supergiant stars. 


(a) Hot dwarfs 

0.13 < 0 < 0.69 

2.75 < log y < 4.37 

polynomial fit 


II 

CO 

1.010 ± 0.333 

ctyp = 0.193 

02 

-40.80 ± 5.24 

Cstd = 0.641 

03 

67.12 ± 7.29 

= 0.85 

(b) Hot giants 

0.14 < 0 < 0.69 

0.69 < log g < 3.01 

polynomial fit 


A = 26 

CO 

-0.4492 ± 0.2930 

ctyp = 0.148 

03 

20.46 ± 1.64 

Cstd = 0.612 



= 0.93 

(c) Warm stars 

0.50 < 0 < 0.90 

0.40 < log g < 4.53 

polynomial fit 


N = 193 

CO 

-27.12 ± 5.62 

o-typ = 0.192 

log y 

18.27 ± 5.69 

Cstd = 0.649 

0 [Fe/H] 

36.08 ± 6.16 

= 0.94 

0 log g 

-29.44 ± 7.25 


02 

127.8 ± 17.2 


log^y 

-2.932 ± 0.983 


[Fe/H] 2 

3.611 ± 1.256 


log g [Fe/H] 

-5.038 ± 1.051 


03 

-73.18 ± 10.86 


logoff [Fe/H] 

0.4446 ± 0.1701 


02 [Fe/H] 

-20.51 ± 4.69 


0 log2g 

4.406 ± 1.239 


log g [Fe/H] 2 

-0.8168 ± 0.3179 


(c’) Warm giants 

0.13 < 0 < 1.10 

0.10 < log g < 3.10 

metal-poor 


[Fe/H] <-0.25 

polynomial fit 


N = 85 

CO 

-0.1213 ± 2.277 

crtyp = 0.230 

0 

8.780 ± 2.056 

Cstd = 0.637 

log y 

-0.3118 ± 0.4345 

= 0.89 

0 [Fe/H] 

2.396 ± 0.337 


(d) Cool stars 

0.70 < 0 < 1.30 

0.00 < log g < 4.85 

polynomial fit 


N = 551 

CO 

-70.87 ± 14.56 

Ctyp = 0.191 

0 

302.2 ± 44.0 

Cstd = 0.540 

log y 

-20.58 ± 2.11 

= 0.95 

[Fe/H] 

57.60 ± 11.20 


0 [Fe/H] 

-80.81 ± 20.67 


0 log g 

12.77 ± 1.88 


02 

-312.5 ± 44.2 


log^y 

3.514 ± 0.453 


[Fe/H] 3 

2.314 ± 0.888 


log g [Fe/H] 

-5.789 ± 0.971 


03 

99.75 ± 14.51 


log^y 

-0.1520 ± 0.0323 


log2g [Fe/H] 

0.1315 ± 0.0738 


02 [Fe/H] 

29.58 ± 9.29 


0 [Fe/H] 2 

-2.103 ± 0.916 


0 log2g 

-1.674 ± 0.329 


0 log g [Fe/H] 

4.071 ± 0.886 



Table 4 - continued 


(e) Cool giants 

1.00 < 0 < 1.40 0.00 < log g < 3.50 

polynomial fit 

N = 287 

CO 

368.9 ± 146.4 atyp = 0.182 

0 

-884.9 ± 382.9 o-^td = 0.529 

log y 

-7.260 ± 2.101 r2 = 0.91 

[Fe/H] 

10.40 ± 2.96 

0 [Fe/H] 

-5.575 ± 2.477 

02 

742.7 ± 333.4 

log^y 

0.6285 ± 0.1689 

log g [Fe/H] 

-0.8552 ± 0.2305 

02 log g 

2.389 ± 1.312 

03 

-209.8 ± 96.4 

(f) Cold dwarfs 

1.07 < 0 < 1.84 4.45 < log g < 5.13 

polynomial fit 

N = 23 

CO 

-104.7 ± 23.8 o-typ = 0.200 

0 

245.1 ± 51.1 o-std = 0.290 

02 

-173.2 ± 36.0 r2 = 0.99 

03 

38.71 ± 8.35 

(g) Cold giants 

1.30 < 0 < 1.73 0.00 < log g < 1.60 

exponential fit 

cte. =2.0 N = 27 

CO 

-29.66 ± 14.28 crtyp = 0.127 

0 

48.10 ± 19.76 o-atd = 0.888 

02 

-18.21 ± 6.80 = 0.83 


0.90 for -1.25 < [Fe/H] < -0.75, and 0.65 <6 < 0.85 for 
—0.75 < [Fe/H] < —0.25). Therefore, Fig. ^ only illustrates 
one particular case (dashed line). 

It is important to note that, while the lines in Figure ^ 
correspond to projections of the fitting functions at par¬ 
ticular values of logg and [Fe/H], the plotted stars span a 
range of atmospheric parameters around these central val¬ 
ues. Therefore, it is not expected that the lines ht all the 
points in the plots. However, there are still some curves 
which do not appear to be well constrained by the observa¬ 
tions in certain regions of the parameter space. For instance, 
there are no stars for 6 < 0.5 in Figure ^ (logg = 1.0) and 
for 6 > 1.3 in Figure ^b (logy = 2.6). Stars with these pa¬ 
rameters do not exist (see Fig. P), and therefore they will 
never be required by the stellar Mpulations models. The 
curves at [Fe/H]= -1-0.5 (Figure]^) and at logy = 0.0 
(Figure ii) for 0.50 < 9 < 1.0 deserve a word of cau¬ 
tion. The problem here is that high abundance, relatively 
hot supergiants are indeed rare and, therefore, uncertainties 
in the CaT* predictions for this region may be present. Note 
that, since the htting functions have been derived using the 
whole sample of stars within each parameter box, the CaT* 
predictions for high metallicity, hot supergiants are partly 
driven by extrapolations of the functional dependence deter¬ 
mined from stars of lower abundances and higher gravities. 
The contribution of this kind of stars (hot supergiants with 
metallicities above solar) at the near-IR spectra of old galax¬ 
ies is not important, although it would not be the case for 
very young stellar populations. 
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Table 5. Coefficients and statistical data of the local fitting func¬ 
tions for the index PaT in each range of atmospheric parameters. 
The term “giant” also includes supergiant stars. 


(a) Hot dwarfs 

0.13 < 0 < 0.69 

2.75 < log g < 4.37 

polynomial fit 


A = 46 

CO 

-6.415 ± 0.792 

Ctyp = 0.165 

e 

54.71 ± 3.39 

Cstd = 0.911 

02 log g 

-4.283 ± 1.091 

= 0.91 

03 

-55.44 ± 5.66 


(b) Hot giants 

0.14 < 0 < 0.69 

0.69 < log g < 3.01 

polynomial fit 


N = 29 

CO 

1.601 ± 1.512 

crtyp = 0.123 

02 

44.81 ± 21.41 

(Zstd = 1.310 

03 

-46.88 ± 29.20 

r2 = 0.64 

(c) Warm stars 

0.50 < 0 < 0.90 

0.40 < log g < 4.53 

polynomial fit 


N = 193 

CO 

-149.0 ± 21.7 

(Ttyp = 0.170 

0 

651.2 ± 88.4 

(Zstd = 0.553 

log g 

18.74 ± 3.83 

= 0.96 

0 [Fe/H] 

-4.651 ± 1.766 


0 log g 

-56.40 ± 10.92 


02 

-840.9 ± 124.0 


02 log g 

39.26 ± 7.62 


03 

338.0 ± 60.0 


02 [Fe/H] 

5.880 ± 2.139 


(d) Cool stars 

0.70 < 0 < 1.30 

0.00 < log g < 4.85 

polynomial fit 


N = 551 

CO 

177.1 ± 12.0 

(7typ = 0.171 

0 

-445.5 ± 33.8 

Cstd = 0.296 

log g 

-16.00 ± 1.38 

= 0.89 

0 log g 

29.20 ± 2.82 


02 

372.9 ± 31.2 


[Fe/H] = 

-0.3610 ± 0.1033 


02 log g 

-13.33 ± 1.43 


03 

-103.3 ± 9.50 


[Fe/H]3 

-0.1080 ± 0.0453 


(e) Very cool stars 

1.00 < 0 < 1.30 

0.00 < log g < 4.85 

polynomial fit 


N = 300 

CO 

0.7600 ± 0.0343 

Ctyp = 0.166 

[Fe/H] 

0.1743 ± 0.0550 

(Zstd = 0.251 



= 0.13 

(f) Cold dwarfs 

1.07 < 0 < 1.84 

4.45 < log g < 5.13 

polynomial fit 


N = 23 

CO 

1.200 ± 0.208 

(7typ = 0.177 

03 

-0.3094 ± 0.0679 

Cstd = 0.374 



= 0.63 

(g) Cold giants 

1.25 < 0 < 1.73 

0.00 < log g < 1.65 

polynomial fit 


II 

CO 

344.3 ± 183.8 

crtyp = 0.112 

0 

-740.2 ± 374.5 

cTstd = 0.351 

02 

527.0 ± 253.0 

= 0.79 

03 

-123.8 ± 56.7 




Figure 8. PaT values for all the stars around solar metallicity 
(—0.25 < [Fe/H] < -1-0.25). The curves correspond to general fit¬ 
ting functions computed for [Fe/H] = 0.0 and different values of 
gravity (logg = 0.0, 1.0, 2.0, 3.0, 4.0 and 5.0, from the thickest 
to the thinnest line). Codes and relative sizes of the star sym¬ 
bols (indicating, respectively, metallicity and gravity ranges) are 
explained in Fig. ^a. 

3.2.2 PaT fitting functions 

As in the case of the CaT*, the PaT index in hot and cold 
stars shows a different trend for dwarfs and giants. Once 
again, we designed proper boxes for each group (a, b, f and 
g in the lower panel of Fig. fe) and derived the local fitting 
functions given in Tables ^f and Figure ^ shows 

the general htting functions for solar metallicity and differ¬ 
ent gravities. For hot stars, PaT increases as temperature 
decreases, reaching a maximum at a temperature that de¬ 
pends on the luminosity class. To reproduce the dependence 
on gravity of the hot dwarfs, a term in log g was deamed 
necessary. At the cold temperature end, the PaT in giant 
stars exhibits a bump due to the presence of strong TiO 
absorptions falling into the index bandpasses. 

Warm stars follow a decreasing trend with increasing 9 
up to the cool stars which, as expected, attain values around 
zero. Boxes c and d were mainly designed to reproduce the 
maximum and the above trends. For cool stars, we have 
found a mild metallicity dependence in the sense that the 
larger the metallicity the stronger the PaT. Finally, box e 
allows smooth interpolations with the coldest stars. 

3.3 Residuals and error analysis 

We dehne the residual of any star as the difference between 
the observed index and the one predicted by the fitting func¬ 
tions (AT = Tobs — Tpred)- Figure ^ shows the residuals of 
the indices CaT*, PaT and CaT for the whole stellar li¬ 
brary as a function of 9. The residuals do not exhibit sys¬ 
tematic trends with any of the three atmospheric parame¬ 
ters. Star clusters have also been analyzed separately and, 
except for the globular cluster M71, no systematic effects 
have been found. Compared to the predictions of the ht¬ 
ting functions, this cluster shows a signihcant mean offset 
of ACaT* = —0.35 and ACaT = —0.51. The most plausible 
source of these systematic offsets is an error in the atmo¬ 
spheric parameters. A systematic error in the controversial 
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Figure 9. Residuals of the derived fitting functions (Al = Xobs — 
2^pred) against 9. See Fig. for symbol codes. 


reddening of this cluster would affect the adopted values. 
However, given the Call triplet dependence on atmospheric 
parameters in this region of the parameters space, only a 
change in the adopted metallicity could account for the de¬ 
rived offsets. The controversy about the metallicity of M71 
has already been discussed in Section 6.1 of Paper II. In fact, 
these residuals were used as one of the arguments in favour 
of the Carretta & Gratton (1997) metallicity scale (note that 
the value of [Fe/H] = —0.56 established by the Zinn & West 
(1984) scale leads to even larger residuals). According to the 
predictions of the new fitting functions, M71 should have a 
metallicity of [Fe/H]= —0.84 ± 0.06, instead of the adopted 
value of [Fe/H] = —0.70. Nevertheless, we have decided not 
to change the input metallicity of this cluster, keeping the 
Carreta & Gratton scale for all the globular clusters. Note 
that the GaT metallicity dependence is constrained by us- 


Table 6. Statistical data for the general fitting functions of the 
indices CaT*, PaT and CaT. (Tstd- unbiased residual standard 
deviation; crtyp: typical index error for the stars used in the fits; 
r^: determination coefficient. 



^std 

(Jtyp 

r2 

CaT* 

0.534 

0.184 

0.965 

PaT 

0.427 

0.163 

0.965 

CaT 

0.629 

0.197 

0.939 


ing more than 500 field and cluster stars, and a change in 
the adopted metallicity of M71 would not significantly al¬ 
ter the results. An anomalous [Ga/Fe] abundance ratio for 
M71 could also explain the derived offsets. However, using 
the tabulated data from Thevenin (1998), we derive a mean 
[Ga/Fe] ratio for this cluster of 0.28 ± 0.07, which is signifi¬ 
cantly above the expected mean [Ga/Fe] value (0.13) for the 
cluster metallicity (see Figure |lC| ) and, therefore, can not 
explain negative residuals in the indices. 

In order to explore in more detail the reliability of the 
present fitting functions, in Table ^ we list the unbiased 
residual standard deviation from the fits, (Jstd, the typical 
error in the measured indices, crtyp, and the determination 
coefficient, r^, for the 664 stars employed in the computation 
of the general fitting functions (note that, apart from the 23 
stars listed in Table ^ 19 stars with unknown [Fe/H] could 
not be included in the metallicity-dependent fits). It is clear 
that, for the three indices, the computed (Jstd is larger than 
what it should be expected uniquely from index uncertain¬ 
ties ^ee also the partial values of Ustd and crtyp in Tables ^ 
and H) . Since we are confident that the latter errors are reli¬ 
able (see full details on index errors computation in Section 5 
of Paper I), the residuals must be dominated by a different 
error source. In fact, such an additional scatter could be due 
to uncertainties in the input atmospheric parameters of the 
library stars. In order to check this, we have computed how 
the parameters errors translate into uncertainties in the pre¬ 
dicted GaT*. This depends both on the local functional form 
of the fitting function (e.g., a weak dependence on temper¬ 
ature leads to small index errors due to Tgft uncertainties) 
and on the atmospheric parameters range (e.g., hot stars 
have Tefi uncertainties larger than cooler stars). Thus, we 
have carried out the following analysis for different groups 
of stars (listed in Table ^ . 

For each star of the sample we have derived three GaT* 
errors, arising from the corresponding uncertainties in Teff, 
logg and [Fe/H]. As input atmospheric parameters uncer¬ 
tainties we have made use of the values presented in Table 7 
of Paper H. Apart from those, we have used errors of 75 K, 
0.40 dex and 0.15 dex for the effective temperatures, grav¬ 
ities and metallicities taken from Soubiran, Katz & Gayrel 
(1998) (stars coded SKC in Table 6 of Paper H), and 75 K, 
0.05 dex and 0.20 dex for the cluster stars. For each star 
group we have computed a mean GaT* error as a result of 
the uncertainty of each parameter ((JTeff, <Jiog g and (J[Fe/ H]) 
by using the input parameters errors for all the individual 
stars. Finally, an estimate of the total expected error due to 
atmospheric parameters ((Jpar) is computed as the quadratic 
sum of the three previous errors. 

Table ^ presents a comparison of these expected er¬ 
rors with the unexplained residual errors of the fits (ures = 
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Table 7. Uncertainties of the CaT* fitting functions for different groups of stars, and mean CaT* errors due to uncertainties in the input 
atmospheric parameters. N: number of stars. cTsta: unbiased residual standard deviation of the fit. crtypi typical observational CaT* error 
for the subsample of stars. Cres : residual error (derived as ^std ~ ^typ) * 9 ^[Fe/H] • niean CaT* errors due to uncertainties 

in the input Teff, \ogg and [Fe/H]. cTpar^ total error due to atmospheric parameters (quadratic sum of the three previous errors). 



N 

^std 

Ctyp 

*^res 


aiogg 

'^[Fe/H] 

^par 

Open clusters (Coma, Hyades, M67, NGC188, NGC7789) 

93 

0.54 

0.28 

0.46 

0.24 

0.07 

0.42 

0.49 

Globular clusters (M3, M5, MIO, M13, M71, M92, NGC6171) 

52 

0.84 

0.76 

0.36 

0.08 

0.11 

0.54 

0.56 

Field dwarfs 

243 

0.45 

0.18 

0.41 

0.23 

0.15 

0.23 

0.36 

Field giants 

201 

0.46 

0.16 

0.43 

0.15 

0.26 

0.23 

0.38 

Field supergiants 

75 

0.84 

0.17 

0.82 

0.47 

0.46 

0.52 

0.84 

Hot stars (0.13 < 6 < 0.69) 

68 

0.66 

0.17 

0.64 

0.81 

0.23 

0.42 

0.94 

Intermediate stars (0.69 < 0 < 1.30) 

555 

0.48 

0.19 

0.45 

0.12 

0.19 

0.31 

0.38 

Cold stars (1.30 < 9 < 1.84) 

41 

0.71 

0.14 

0.70 

0.69 

0.20 

0.21 

0.75 

All 

664 

0.53 

0.18 

0.50 

0.22 

0.20 

0.31 

0.43 



Figure 10. [Ca/Fe] ratios against [Fe/H] for 216 stars of the li¬ 
brary. Different symbols and sizes indicate metallicities and grav¬ 
ities as in Fig. The relationship [Ca/Fejg = 0.013(di0.012) — 
0.173(ih0.011) [Fe/H] (solid line) is a least-squares fit to the data 
and confirms an intrinsic anticorrelation between both abundance 
ratios. 

It is clear from this table that, in most cases, 
the additional scatter (ares) is comparable to apar, that is, it 
can be explained from uncertainties in the input atmospheric 
parameters. Using an F test of comparison of variances we 
have checked that Ures is not significantly larger than Upar, 
except for the intermediate stars. This means that, a mi¬ 
nor, additional error source may still be needed to account 
for the observed residuals of field dwarfs and giants with 
intermediate temperatures. 

Finally, since the aim of this paper is to predict reli¬ 
able index values for any combination of input atmospheric 
parameters, we have also computed, making use of the co- 
variance matrices of the fits, the random errors in such pre¬ 
dictions. These uncertainties are given in Table ^ for some 
representative values of input parameters. Note that, as is 
expected, the absolute errors are larger for cold giants and 
supergiants and, in general, increase as the metallicity de¬ 
parts from the solar value. 

3.4 Ca/Fe abundance ratios 

Since CaT* is a pure Can index, the parameter governing 
the metallicity dependence should be the actual [Ca/H] ra- 
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Figure 11. Residuals of the CaT* fitting functions (observed 
— predicted) versus residuals of the linear [Ca/Fel relationship. 
Symbols and sizes for the 216 stars are as in Fig. ph. 


tio, rather than the classical [Fe/H]. Since, in spite of this, 
we have been using [Fe/H] as our metallicity parameter, we 
have investigated whether the residuals of the CaT* fitting 
functions are correlated with the individual Ca/Fe abun¬ 
dance ratios. Note that, if that were the case, this could 
account for part of the unexplained residuals reported in 
the previous section. Using a subsample of 216 library stars, 
for which we have compiled [Ca/Fe] data from the literature 
(Gratton & Sneden 1987; Pilachowski, Sneden & Kraft 1996; 
Nissen & Schuster 1997; Thevenin 1998), we have not found 
any evidence for such an expected correlation. This agrees 
with the results of ITD, which, after including a [Ca/Fe] 
term in their fitting functions, concluded that this relative 
abundance did not practically affect the value of CaT. In 
fact, a proper error analysis demonstrates that their derived 
[Ca/Fe] coefficient is not statistically different from zero. 

To investigate this issue in more detail, in Figure 
we plot the [Ca/Fe] relative abundances versus the assumed 
metallicities in [Fe/H]. As it has been previously reported 
(e.g. Edvardsson et al. 1993), there exists a clear anticor¬ 
relation between both abundance ratios for Galactic stars. 
A least-squares linear fit to all the stars in this plot gives 
[Ca/Fe]o = 0.013(±0.012)-0.173(±0.011) [Fe/H]. The exis¬ 
tence of this relation indicates that the actual calcium abun¬ 
dance has been implicitly taken into account in the fitting 
functions through the adopted [Fe/H] and, therefore, a sys- 
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Table 8. Absolute errors in the fitting functions predictions for different values of the atmospheric parameters. Input logg values varying 
with effective temperature for dwarfs, giants and supergiants have been taken from Lang (1991). Since, for extreme temperatures, the 
fitting functions do not depend on metallicity, no [Fe/H] value has been adopted for 15000 K and 3200 K. This is also the case for some 
values at 3500 K. 





dwarfs 



giants 


supergiants 

Teff 

[Fe/H] 

ACaT* 

APaT ACaT 

ACaT* APaT ACaT 

ACaT* APaT ACaT 

15000 


0.17 

0.25 

0.29 

0.17 

0.25 

0.29 

0.24 

0.60 

0.61 

8000 

+0.5 

0.39 

0.29 

0.47 

0.45 

0.25 

0.51 

0.60 

0.32 

0.67 

8000 

0.0 

0.28 

0.24 

0.36 

0.30 

0.20 

0.36 

0.29 

0.28 

0.39 

8000 

-1.0 

0.40 

0.33 

0.50 

0.37 

0.32 

0.47 

0.80 

0.39 

0.87 

8000 

-2.0 

0.76 

0.54 

0.91 

0.77 

0.54 

0.92 

0.81 

0.59 

0.98 

6000 

+0.5 

0.25 

0.09 

0.26 

0.38 

0.11 

0.39 

0.60 

0.20 

0.63 

6000 

0.0 

0.10 

0.06 

0.11 

0.22 

0.09 

0.23 

0.37 

0.20 

0.41 

6000 

-1.0 

0.17 

0.09 

0.19 

0.34 

0.11 

0.35 

0.50 

0.21 

0.54 

6000 

-2.0 

0.37 

0.16 

0.40 

0.58 

0.18 

0.60 

0.68 

0.25 

0.72 

5000 

+0.5 

0.28 

0.07 

0.29 

0.19 

0.04 

0.20 

0.31 

0.08 

0.32 

5000 

0.0 

0.17 

0.07 

0.18 

0.10 

0.04 

0.11 

0.18 

0.08 

0.20 

5000 

-1.0 

0.19 

0.08 

0.20 

0.14 

0.06 

0.15 

0.28 

0.09 

0.29 

5000 

-2.0 

0.39 

0.12 

0.40 

0.29 

0.10 

0.31 

0.41 

0.11 

0.42 

4000 

+0.5 

0.11 

0.08 

0.13 

0.41 

0.13 

0.42 

0.66 

0.13 

0.68 

4000 

0.0 

0.11 

0.07 

0.13 

0.26 

0.12 

0.28 

0.54 

0.12 

0.55 

4000 

-1.0 

0.11 

0.07 

0.13 

0.38 

0.13 

0.40 

0.40 

0.13 

0.42 

4000 

-2.0 

0.11 

0.10 

0.15 

0.74 

0.16 

0.75 

0.54 

0.16 

0.56 

3500 

+0.5 

0.13 

0.10 

0.16 

0.97 

0.19 

0.98 

0.99 

0.19 

1.00 

3500 

0.0 

0.13 

0.10 

0.16 

1.00 

0.19 

1.02 

1.00 

0.19 

1.01 

3500 

-1.0 

0.13 

0.10 

0.16 

1.13 

0.19 

1.14 

1.07 

0.19 

1.09 

3500 

-2.0 

0.13 

0.10 

0.16 

1.30 

0.19 

1.31 

1.20 

0.19 

1.22 

3200 


0.19 

0.13 

0.23 

0.37 

0.22 

0.42 

0.37 

0.22 

0.42 


tematic trend between the fitting functions residuals and 
the Ca/Fe ratios should not be expected. However we could 
still have a correlation between the residuals from the above 
relation (A[Ca/Fe] = [Ca/Fe] — [Ca/Fe]Q), and the fitting 
functions residuals. Figure]^ shows that these residuals are 
not correlated at all. Thus, we conclude that the new fitting 
functions do not need any additional term to account for 
Ca overabundances within the library stars. Note however, 
that the fitting functions in their present form (i.e. using 
[Fe/H] as the metallicity indicator) implicitly include the 
chemical enrichment of the solar neighborhood, although, 
in principle, they could be corrected, by introducing a new 
[Ca/Fe]-[Fe/H] relation, to give predictions for other enrich¬ 
ment scenarios. 


4 COMPARISON WITH PREVIOUS FITTING 
FUNCTIONS 

In this section, we present a comparative analysis of the 
prev ious fitting functions by DTT, ITD and JCJ (see Sec- 
and Table & and those presented in this paper. DTT 


tion 2.2 


and ITD made use of empirical stellar libraries of 106 and 
67 stars repectively, and derived fitting functions for their 
calcium indices. On the other hand, using NLTE models, 
JCJ derived fitting functions for the theoretical equivalent 
widths of the Ca ii lines. 

In order to avoid extrapolations of the previous fit¬ 
ting functions, the range of effective temperatures in which 


comparisons have been performed was constrained to the 
one approximately spanned by the three previous works 
(0.75 < 0 < 1.25). Also, since the empirical papers by 
DTT and ITD do not correct for the Paschen contamina¬ 
tion, their predictions are compared with our CaT fitting 
functions. For the comparison with JCJ we use the CaT* 
fitting functions to compare with their theoretical Ca ll pre¬ 
dictions, since they are based on model atmospheres which 
do not include the contamination by other elements. 

Prior to any comparison, the previous fitting functions 
have been corrected from systematic offsets arising from dif¬ 
ferences in the index definitions. For the predictions by DTT 
and ITD, we have applied the calibrations given in Table 8 of 
Paper I, which were specially designed to convert the values 
of the literature indices to our system. Although JCJ do not 
define line-strength indices, they provide fitting functions 
for the total true equivalent width of the two strongest Ca ll 
lines. We have increased their predicted equivalent widths by 
21.1(±0.3) percent. Such a fraction, which has been empir¬ 
ically derived using all the library stars within the range of 
effective temperatures given above, quantifies the strength 
of the weakest line of the triplet relative to the sum of the 
two other ones. Hence, after these corrections, differences 
between the predictions of the fitting functions will only 
arise from differences in the input data quality (in the case 
of the empirical papers, this will include the observational 
errors in the indices, the uncertainties in the atmospheric 
parameters and the number of calibrating stars), the range 
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of atmospheric parameters spanned by the library, and the 
mathematical procedure to compute the fits. 

Figure shows the present Call fitting functions 
(dashed lines) together with those derived by DTT, ITD and 
JCJ (solid lines). The predicted gravity effects are shown in 
the upper panels (a, b, c), whereas mid (d, e, f) and lower 
panels (g, h, i) show, respectively, the predicted dependences 
on metallicity for dwarf and giants. 

DTT found a biparametrical, linear Call behaviour 
with metallicity and gravity, but no dep ende nce on tem¬ 
perature (see also the discussion in Section 2.2). Concerning 
the gravity effect, DTT predict a milder dependence than 
what our data shows, as it is apparent from Fig. [l2^ . Also, 
there are important differences at the other temperature 
ranges, although metallicity effects agree with our predic¬ 
tions, especially for warm dwarfs (Figs. p^) and cool gi¬ 
ants (Figs. [l^). We must note that most differences be¬ 
tween DTT and our predictions arise from: i) the lack of 
second order and loggf-[Fe/H] cross-terms in their derived 
functions (a reanalysis of their data shows that these terms 
were indeed statistically significant), and ii) the absence of 
an effective temperature dependence. This last discrepancy 
is mainly due the lack of cold {9 > 1.05) dwarfs in their 
sample (see the effect in panel d) and to the contamination 
by Paschen lines (panel g). 

The fitting functions by ITD include terms in effec¬ 
tive temperature and several cross-terms in the parameters. 
They concluded that the gravity effect is not as important as 
the dependence on metallicity. Fig. |l^ shows that their lin¬ 
ear, extremely weak, gravity dependence, does not agre e at 
all with our conclusions. As it was discussed in Section |2.2[ 
this is mostly due to the scarcity of supergiants in their 
sample and to the fact that most of them have very low 
metallicities. On the other hand, the metallicity dependence 
for dwarfs and cool giants roughly matches with our pre¬ 
diction^ although it is again very different for warm giants 
(Figs. |i|e,h). 

JCJ found a complex dependence on the three atmo¬ 
spheric parameters. Unlike the two previous empirical works, 
there is a good agreement between the qualitative behaviour 
of the parametrized gravity dependence (Fig. [l^), in the 
sense that the lower the gravity, the larger the gravity effect. 
Also, they find a metallicity dependence with a functional 
form similar to our predictions, although their values are 
more sensitive to metallicity than ours, especially for dwarfs 
(Figs.@f,i). 

We conclude that, in the temperature range of compar¬ 
isons, the gravity dependence of the Call triplet predicted 
by JCJ, rather than that of the empirical papers, roughly 
agrees with our results. On the other hand, for certain tem¬ 
perature ranges, DTT and ITD better reproduce the metal¬ 
licity dependence found in this work. 


5 CONCLUSIONS 

In this paper, we have analysed the behaviour of the Call 
triplet strength in the spectra of stars of different spectral 
types and luminosities, by means of a new near-IR stellar 
library (presented in Paper I) with a wide coverage of atmo¬ 
spheric parameters (Paper II). We have derived empirical 
fitting functions which can be easily implemented into stel- 
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lar populations models. Readers interested in including these 
functions into their population synthesis codes can make use 
of the FORTRAN subroutine referred to in Section |^. 

Other results of this work are the following: i) We find a 
complex behaviour of the Ca ll strength as a function of the 
three main atmospheric parameters (Teg, logg and [Fe/H]). 
For hot and cold stars, effective temperature and luminosity 
class are the main driving parameters, whereas, in the mid¬ 
temperature regime, all three parameters play an important 
role (see Fig. M). ii) The residuals of the fitting functions 
arise mainly from uncertainties in the input atmospheric 
parameters, iii) We do not find any correlation between these 
residuals and [Ca/Fe] abundance ratios, iv) A comparison 
with the fitting functions derived in the literature reveals 
striking differences in the predicted gravity and metallicity 
dependences for some regions of the parameter space (see 
Fig. 1^. 

Basically, compared to previous fitting functions for the 
Call triplet, the main advantage of our predictions is that 
they include the whole range of effective temperatures. This 
is extremely important since the Ca ll behaviour in the tem¬ 
perature ranges explored by the previous works can not be 
extrapolated at all to hot (6 ^ 0.75) or cold stars (9 ^ 1.35), 
as it is readily seen in, for instance. Fig. ^ We should not 
forget that these cold stars constitute a significant contribu¬ 
tion to the integrated near-IR spectra of most stellar popu¬ 
lations (to old populations through the role of cold dwarfs 
and giants, and to relatively young populations through the 
expected contribution of AGB stars; see Paper IV). This 
must indeed be remarked since, if the fitting functions by the 
previous works are extrapolated, stellar populations synthe¬ 
sis models will derive unreliable integrated Call strengths 
which can not be compared with observations. 

In the case of synthesis models that make use of the 
theoretical fitting functions of JCJ, the comparison of the 
predictions with the observed Call triplet in galactic spec¬ 
tra is even more uncertain since: i) JCJ predictions do not 
include the contamination of Paschen lines, which affects 
the measurement of the observed index if the contribution 
of warm stars is not negligible; and ii) JCJ provides true 
equivalent widths, whereas line-strength indices in systems 
like that of DTT are used for the observational data. Al¬ 
though JCJ measure the strengths using the central band- 
passes of DTT, these measurements can not be compared 
to indices in this system, since the local pseudo-continuum 
(traced between the sidebands) is far from being a true con¬ 
tinuum (especially for hot and cold stars; see Paper I). 

As a final conclusion, we want to emphasize the impor¬ 
tance of using an empirical library with a wide range of ac¬ 
curate atmospheric parameters and a sufficiently large num¬ 
ber of calibrating stars for deriving reliable empirical fitting 
functions. Also, well defined index definitions and objective 
fitting procedures are a critical factor. To summarize, the 
functions presented in this paper alleviate these problems 
and should lead to more reliable comparisons between mod¬ 
els and observations than in the past. The first predictions 
of population synthesis models using these results will be 
presented in Paper IV. 
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Figure 12. Comparison with previous fitting functions. Panels a, b and c display, respectively, the predictions by DTT, ITD and JCJ 
for solar metallicity stars and different gravities (log^ = 0.5, 1.5, 2.5, 3.5, and 4.5). The thicker the line, the lower the gravity. Solid lines 
correspond to the fitting functions by the above authors, whilst our predictions are given in dashed lines. Analogously, the predicted 
dependence on metallicity is also compared for dwarfs (log^ = 4.5), in panels d, e and f, and giants (log^ = 2.5), in panels g, h and i. 
The different lines represent, from top to bottom, the metallicities [Fe/H] = +0.5, 0.0, -0.5, —1.0, —1.5, and -2.0, for DTT and ITD, and 
[Fe/H] = +0.2, 0.0, -0.5, and —1.0 for JCJ. Note that, in all cases, the predictions by DTT and ITD are compared with the CaT fitting 
functions, whereas the theoretical work by JCJ is compared with CaT*. 
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